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Abstract 



In this paper, for an operator defined by the action of an M-th order differ- 
ential operator with rational-type coefficients on the function space L^j^^^(M) := 

{/ : measurable] ||/||(fco) < 00} with norm H/H^^^) := J \f{x)\'^{x'^ + l)^''dx 
(feo S Z), we prove the regularity (continuity and differentiability up to M 
times) of the eigenfunctions of its closure (with respect to the graph norm), 
except at singular points of the corresponding ordinary differential equation 
without any assumptions for the Sobolev space, i.e., without any assumptions 
about the m-th order derivatives of the eigenfunctions with m = 1, 2, . . . M—1. 
(For the special case of ko = 0, we prove this regularity for the usual L^(M).) 
Especially, we show a one-to-one correspondence between the eigenfunctions 
of its closure and the solutions in C^(M) H L^^^-|(]R) of the corresponding dif- 
ferential equation under the condition above when there is no singular point 
for this differential equation. This one-to-one correspondence is shown in the 
basic framework of an algorithm proposed in our preceding paper, which can 
determine all solutions in C*^ n L^^^j(M) of the ordinary differential equation 
then. 
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1 Introduction 



When we treat the eigenfunction problem of the closure of an M-th order differen- 
tial operator on a Hilbert space with a certain boundary condition, we should be 
careful to distinguish this problem from the problem of finding solutions in the space 
of M-times differentiable functions C*^(R) of the differential equation described by 
this differential operator definable only in C^'^M.). If a solution to the latter problem 
belongs to the Hilbert space above and satisfies the boundary condition, it is always 
an eigenfunction of the former problem from the definition. However, it is not neces- 
sarily the case that eigenfunctions of the former problem belong to C*^(R). Hence, 
the regularity (continuity and differentiability up to M times) of the eigenfunctions 
of the former problem should be examined carefully. 

In the theory of elliptic operators pQ, this problem has been discussed under as- 
sumptions for the Sobolev space, i.e., the assumption that the m-th order derivatives 
of the eigenfunction with m = 1, 2, ... M — 1 belong to L^-space. These assump- 
tions are often required for the validity of numerical methods that solve differential 
equations by projection to finite dimensional subspaces (Ritz-Galerkin and Petrov- 
Galerkin methods [2] [3], for example). 

On the other hand, in this paper, for a class of Hilbert spaces containing L^(M), 
we will discuss the regularity problem above under several conditions, without any 
assumptions concerning the m-th order derivatives of the eigenfunction. The condi- 
tion in our discussion is that the differential operator has rational coefficient functions 
and its characteristic equation (eigenvalue equation) has no singular point. Under 
this condition, we prove that the regularity above is always guaranteed. 

This discussion can be generalized, even for a differential operator (with rational 
coefficient functions) whose characteristic equation has singular points, by excluding 
only the regularity at the singular points of the characteristic equation. Especially 
for Fuchsian-type differential operators, we give a stronger statement than general 
cases. 

The proof is based on a one-to-one correspondence between the 'regular' solutions 
in the Hilbert space of the differential equation and the square-summable number- 
sequence solutions of simultaneous linear equations described by a kind of matrix 
representation of the action of the differential operator, which is guaranteed under 
several conditions. In this paper, we will clarify how we can show regularity using 
this one-to-one correspondence. 

The contents of this paper are as follows: Section [2] introduces the basic frame- 
work used for the proof. Firstly, in Subsection 12.11 we clarify precisely what has to 
be proved. Next, Subsection 12.21 provides a more general framework in which the 
regularity problem can be discussed, and it shows the conditions that are required 
for the base of this framework. In Subsection 12. 3[ In order to treat the argument 
given in Subsection 12. we give a concrete structure for the general theory given in 
Subsection 12.21 Subsection 12.41 explain how to apply general theory given in Subsec- 
tion 12.21 to the concrete structure given in Subsection 12.31 for showing the statement 
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given in Subsection 12.11 Section [3] is devoted to a proof of a property mentioned in 
Subsection 12.31 Section H] is devoted to a proof of a theorem mentioned in Subsection 

m 

2 Basic framework of this paper 

2.1 'Regularity' of eigenfunctions to be shown 

In this subsection, we rigorously describe the regularity problem to be solved in this 
paper. In this paper, we treat the differential operator 

M 
m=0 

on the space of M-times differentiable functions C^^(]R). 

In order to treat the ODE R{x,-^)f{x) = A/(x) using functional analysis, we 
have to define the differential operator in a complete function space. 

In the present paper, we focus on the function space L^^^-|(]R), which is defined 

by 

L2,^)(M) := {/ : measurable | < oo} (2) 

/oo poo 
and norm ll/lljfc^,) = / \f{x)\^{x^+ 
oo J — oo 

l)'"^dx, parametrized by ko G Z. Here, in the special case of = 0, L^q^(M) is iden- 
tical to the usual L^(M). Then, the operator Ar,l^^^ ^(m) is defined by the action of 
R(x, ^) with domain 

^(^H,L^,^,(R)) := {/ e C'iR) n 4,„)(M) I Rix, ±)f e 4„)(M)}, (3) 

and its closure ^_r,l2^ j(r) with respect to the graph norm [1]. 

In general, an eigenfunction of the closed extension of the given differential oper- 
ator does not necessarily yield a solution of the ODE _R(x, -^)f{x) = Xf{x). This is 
because there is a possibility that the eigenfunction is not an M-times differentiable 
function. This problem is called the regularity problem for a differential operator. 

In the present paper, we prove that the eigenfunction of the operator A^^2 m-, 

always does yield a solution of the ODE R{x, -^)f{x) = A/(x). Since this problem 
depends on the singularity of the differential equation, we need the following def- 
inition. A real number a; G M is called a singular point of a differential operator 
R{x, ^) if X is a singular point of the differential equation R(x, -£:)f{x) = A/(x) 
for some real number A. Indeed, the above definition does not depend on A. When 
the coefficient functions of P{x,^) := Ylim=QPTn{x){-^)"^ are polynomial, the set 
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of its singular points equals the set of zero points of its coefficient function pm of 
the highest degree, which is written as p^/ (0). When the coefficient functions rm{x) 
[m = 0,1,... M) of R{x,-£:) are rational functions, we denote the least common 
multiple of the denominators of rm{x) by l{x). Then, the set of its singular points 
equals the set of zero points of I-tm^x) := /(x)rAf (x), which is written as (/-rAf )~-^(0). 
First, we obtain the following theorem. 

Theorem 2.1 Assume that a differential operator R{x, ^) has no singular points 
and its coefficient functions rm{x) (m = 0, 1, . . . M) are rational functions. Then, an 
element f G L^^^^(M) is an eigenfunction of Aj^ i2^ ^(jg) with the eigenvalue A if and 

only if f belongs to C^^(R) n ^ka)(^) satisfies the ODE R{x, £)f{x) = Xf{x). 

This theorem is one of the main statements to be proved in this paper. When the 
differential operator R{x, ^) has singular points, the following extension holds. 

Theorem 2.2 When the coefficient functions rm{x) (m = 0, 1, . . . M) of a differ- 
ential operator R{x , are rational functions, then any eigenfunction of A^ i^2^ ^(jg) 

belongs to C^^(R \ (/ ■ rM)-\0)) for any integer /co- 
in the next subsection, we will give a more general argument, which includes 
Theorem 12.11 as a special case. 

2.2 Regularity in a more general framework 

In this subsection, we treat the regularity problem in a general Hilbert space H of 
functions on the real line M. That is, we give three conditions equivalent to the 
solution of the ODE P{x,£)f{x) = Xf{x) in a general Hilbert function space H, 
where we convert the ODE to square-summable solutions of a matrix- vector equation 
(simultaneous linear equations) defined in the following general framework. 

Now, we introduce another general Hilbert function space Ti^ as a Hilbert function 
space on M which contains (as a subset) the original Hilbert function space T-t. In 
general, the inner product of "H is distinct from the inner product of "H^, whereas 
"H is a subset of "H^. By treating the differential operator as an operator from "H 
to "H^, we are able to utilize a 'matrix representation' of the ODE with respect 
to appropriate basis systems. The key point of the method that we present is the 
difference between the inner products of the spaces "H and Ti^. 

Define the operator Ap^y^ as the action of P{x, ^) with domain 



and its closure Ap^-^ with respect to the graph norm. Next, we introduce an operator 
from "H to Ti^. Define the operator -Bp^A.w.^o ^is the action of 



D{Ap,n) ■■= {/ e n n I P(x, £)f e n}, 



(4) 



P(A;x,£) :=Pix,£)~XI 



(J : identity op.) 



(5) 
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with domain 

D(^p,,,^,^o) {/ e C^(M) n U I P(A; x, ^)/ e (6) 

and its closure Bp^nn^ with respect to the corresponding graph norm || • 11^^ + 

\\PiX;x,£) .\\^,." ' 

In order to using a band-diagonal structure in the close operator Bpx^n,H'^j 
troduce Conditions C1-C3, Cl"^, C2+, and C2.1-C2.3 for the quintuplet consisting 
of the linear differential operator P(A; x, the Hilbert spaces H and H^, and their 
CONSs {e, and {e^}^=o, which is abbreviated to (P(A; x, {e„ }~ o, {e^}^=o)- 

In what follows, (•, •)^o and (•, denote the inner products of and H respec- 
tively. These conditions are shown to hold in several examples for P{x, ^) later. 

CI For any n, e„ belongs to D{Bpx^'^^n'>)- 

/oo 
f{x)g{x)v{x)dx. 
-oo 

C2 There exists an integer io such that {Bpx,n,HO^n, ^m)HO — when \n — 

m| > if). 

/oo 
f{x)g{x)v^{x)dx. 
-oo 

16" I 

C2.1 sup ^ < oo. 

nez+\{o} n 

C2.2 The basis functions {n e Z+) belong to C^(R) and there exists a first-order 
differential operator N{x, ^) = ni{x)-^ + no{x) satisfying (a) and (b) below: 

(a) : The functions rii and no belong to C^"^(]R) 

(b) : There exist real numbers {n e Z+) such that N{x, ^) = AnC^ for 
any n e Z+, and liminf J— ^ > 0. 

n->-oo n 

C2.3 There exists a function a in C°(M) such that "^n e Z+ and e R, |e^(x)| < 
d{x). 

C3 There exists a linear operator Cpx^-^ no with domain D{Cpx;H,H^) from a dense 
subspace of to such that e D{Cp^x,H,Ho) and {Bp^x,H,Hof^(^m)HO = 
{f,Cp,x,H,Hoet)H for / e D{Bp^x,H,Ho)- 

Our main issue is the correspondence between the following two kinds of so- 
lutions One kind of solutions are the square-summable solutions of the system of 
simultaneous linear equations corresponding to the matrix representation : = 
{Bpx^n,H^^n, em)wo. The other kinds of solutions are the solutions of the ODE 
P{x, ^)f{x) = \f{x) in a general Hilbert function space "H. 
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Due to the condition C3, the basis belongs to the domain of the adjoint 
operator Bp^^ In the following two conditions, M denotes the order of ^). 
With defined in C2, define the solution space as a space of number sequences 

QO 

V:={{fn}^=o\Y.^Un = ^ (mGZ+)}. (7) 

n=0 

With this definition, one of the 'equivalent conditions' mentioned above is {fn\^=o ^ 
V n £^(Z+). As is shown later, the following theorem holds. 

Theorem 2.3 When the quintuplet (P(x, £), ^, {e„ }~=o> {en}^=o) satisfies Cl- 
C3, C2+, and C2.1-C2.3, then the relations (i) =^ (ii) =^ {Hi) =^ (iv) holds con- 
cerning the conditions for / G "H/ 

(i) : f e domAp,^ and Ap^n f = Xf ■ 

(ii) : f e dom Bp x n n<> ^'^^ Bp x hm^ / = 
{ill): f en and{{f,en)n}^=o'(^Vnf{Z+). 

(iv): f e C^^R \S)n'H and'^x eR\S, P{x, £)f{x) = \f{x), S is the set of 
singular points of (P(x, ^). 

We can prove a stronger argument than the above theorem with an additional 
condition. 

In this paper, we define the following two types of Fuchsian operators: 

Definition 2.1 P{x, -£) is called 'Fuchsian of Type I' if its all singular points are 
regular singular points of the ODE P{x, -^)f{x) = 0. 

Definition 2.2 P{x, ^) is called 'Fuchsian of Type IF if its all singular points are 
regular singular points of the ODE P{x, -£)f{x) = A/(x) for any complex number A. 

These definitions do not depend on A. From the definitions, a Fuchsin operator of 
type II is always Fuchsian of Type I. The both definitions are equivalent when the 
coefficient function pm{x) of the highest order has no zero point whose multiplicity is 
greater than M. (However, otherwise, they are not always equivalent. For example, 
the differential operator + x is a Fuchsian operator of Type I but it is not of 
type II.) The definition of Fuchsian operator used in another paper [S] of us is 'of 
type r in this paper. In the cases of Fuchsian operators of type II, when CI"*" and 
C2"*" holds. Theorem 2.7 in [5] guarantees the implication (iv) =^ (ii) in Theorem 
12. 3[ we obtain the following theorem: 

Theorem 2.4 When the quintuplet {P{x, -£),n, {e„ }^^q, , {e2}^=o) satisfies Cl- 
C3, C1+, C2^, and C2.1-C2.3 and the differential operator P{x, ^) is Fuchsian of 
type II, then the following conditions are equivalent for f 

{ii): f G dom5p^;,^^^^o and Bp^x,H,HO / = 

(m).- f en and{{f,en)n}^=o^V n£\Z+). 

{iv): f G C^^{R \S)nn and'^x eR\S, P{x, £)f{x) = Xf{x), S is the set of 
singular points of {P{x, -£). 



6 



When a differential operator P{x, ^) has no singular point, this operator is a 
special case of Fuchsian differential operators. In this special case, since the set S is 
empty, the relation (iv) =^ (i) is trivial. Then, we obtain the following theorem. 

Theorem 2.5 When the qmntuplet {P{x, £), {e„ }~=0' {(^n}n=o) satisfies Cl- 
C3, C2^ , and C2.1-C2.3 and the differential operator P{x, ^) is has no zero points, 
then the following conditions are equivalent for f eH: 

(i) : f e dom Ap^n and Ap^^ f = ■ 

(ii) : f G domBp^x,H,no and Bp^x,n,no / = 
(m).- feu and {(/,e„)^}- q eVni^Z-^). 

{tv): f e C^{R) n n and G R, P{x, £)f{x) = Xf{x). 

Here, we explain the structure of proof of Theorem 12.31 The statements (z) =^ {ii) 
under C1-C3 can be shown from the following lemma. 

Lemma 2.1 // / G dom Ap^-^ and Ap^i f = Xf, then f G domi?p;^^^^^o and 
BpxHM'^ / = 0. 

Proof: The inclusion relation "H C "H*^ in the sense of sets implies also that any 
function sequence converging for the norm || • ||^ converges for the norm || • ||^o. 
Hence, from the definitions, domAp-^ = dom(ylp^^ — A/) C domi?p;^^^^^o. Since 
the equality Apy^f = Xf i.e. {Apy^ — XI) f = implies Bp^^-u^yof = 0, this suffices 
for the proof of this lemma. ■ 
The statements {ii) =^ {Hi) under C1-C3 can be shown by application of The- 
orem 2.2 of [5j to the operator Bpxy yo. The remaining part {Hi) =^ {iv) will be 
shown in Section [3^2] under CI, C2, 02"*" and C2.1-C2.3. Therefore, our remaining 
tasks are summarized as follows. 

Task 1 (Subsection 12. 3p Constructing "H, "H^, and their CONSs satisfying Condi- 
tions C1-C3, C1+, C2+, and C2.1-C2.3 for a differential operator P{x, £) 
with polynomial coefficient functions Pm{x)- Check of C1-C2, C1+, C2^ has 
been done in [5]. 

Task 2 (Section [3]) Checking Condition C3 in the construction of "H, "H^, and their 
CONSs given in Subsection 12.31 

Task 3 (Subsection 12. 4p Showing Theorems 12 . 1 1 and 12.21 using Theorem 12.31 and the 
construction of Ti, Ti^, and their CONSs given in Subsection 12.31 Indeed, 
when the differential operator R{x, ^) mentioned in Theorems 12. II and 12.21 has 
polynomial coefficient functions, the argument of these theorems are immediate 
from Theorem 12.31 and the above construction. However, it is not trivial in the 
non-polynomial case. 

Task 4 (Section H]) Showing the relation {Hi) =^ {iv) mentioned in Theorem 12.31 
under CI, C2, C2+ and C2.1-C2.3. 
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2.3 Function spaces and basis systems satisfying the condi- 
tions 

In this subsection, we treat the case satisfying the following: (1) The differential 
operator P{x,-^) has polynomial coefficient functions. (2) l-i = (L^^^^(M), T-L^ = 
A (R)), (3) fco and /cq satisfy <kQ — sq with 

sq := max (degp™ - m) . (8) 

m 

The purpose of this subsection is giving CONSs of L?, -.(M) and (M) such that 

satisfying Conditions C1-C3, CI"'", 02"*", and C2.1-C2.3 with the above conditions. 
First, we introduce basis systems {e„ | n G Z"*"} and {e^ | n G Z"^} of L^j^.^^(R) and 



with 



:= L-¥J -^(-ir^'^H^J (10) 
1 

^k,n{^) ■■= ■^—^i^] {neZ), (11) 



[x + + i 

where [aj denotes the largest integer not greater than a. It is easy to show that this 
function satisfies the following properties. 



^fc,n e -^(fc)(lK) , ^k,n{x) = tpk,-h-k-l{x) and (V'fc.m , i^k,h){k) =T^5fhn ■ (12) 

Moreover, they satisfy the following lemma: 

Lemma 2.2 ^\J^'^k,h | ^ G z| is an orthonormal basis o/L^^.-|(M). 

The orthonormal property is shown by f|T2|) . though the proof of completeness is 
somewhat complicated. Its proof is given in Appendix A of [5J. This lemma guaran- 
tees CI. 

The indices of functions in | ^ko,h | ?i G are bilaterally expressed, while the 
indices of basis functions in {e^ | n G Z"*"} are unilaterally expressed, and they are 
'matched' to one another by the one-to-one mapping defined by (fTOj) . In order to 
avoid confusion between them, in this paper, the integer indices with double dots 
denote the bilateral ones in Z, in contrast to the unilateral ones (without double 
dots) in Z"*". 

Since the mapping n — )■ Uk.n is one-to-one from Z+ to Z, the basis systems 

{cn I n G Z+} and {e^ | n G Z+} are identical to ^'\J^^ko,n \ n ^ z| and | -y/^'^fc<> ji | ^ ^ 
respectively. Hence, from Lemma 12. 2[ we have 
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Theorem 2.6 {e„ | n G Z+} and {e^ | n E Z+} are orthonormal basis systems forTi 
and Ti^ , respectively. 

The 'matched' number iik^n in ( p^Oj) has the property 



LI J + i (A: : even) 



^fc,n + ^ I = <( (13) 
[^J (A; : odd) 

which is used later. 

As well as satisfying the orthogonality property above, they satisfy other orthogonality- 
like relations (w.r.t. other inner products) given in [6], one of which is related to 
su(l, l)-number-states [7]. When A; > 0, as is explained in the paper [6] in detail, 
4^k,h{x) is an 'almost-sinusoidally' oscillating wavepacket with a spindle-shaped en- 
velope \ipk,h{x)\ = (x^ -|- 1) i~ , and its approximation to a sinusoidal wavepacket 
with a Gaussian envelope holds for sufficiently large k with respect to the L^-norm. 

In the following part of this subsection, we show that the quintuplet (P, L^^^j(]R), {y^ Vfco, 

satisfies Conditions C2, C2+, and C2.1-C2.3. However, our proof for C3 requires 
several pages, and it will be given in Section 13.21 after the introduction of a tool for 
it in Section [3l 

Firstly, C2^ is obvious from the definition of {■,-){k'^)- Moreover, the definition 
of ipk,n{x) results in the following lemma: 

Lemma 2.3 '^k,n e C^{R) n Lf,^^ . 

Since |'?/'fcO^(x)| = (x^ + 1) 2 holds for any real number x, C2.3 is obvious for 

r— -fcj^+i 

d{x) = \ - (x^ + 1) 2 .In order to show C2.2, we focus on the equality: 

- ^ [{x^ + 1) + (A; + l)x^ iJk,n{x) ={n + i^kd^)- (14) 

Then, the operator A^(x, ^) := — f (x^ + 1) ^ + (^0 + 1)^ satisfies the eigen equation 
A^(x, ^) e2(x) = A„e2(x), where A„ := ri^o „ + -^2^- Since ffT5]) implies the inequality 

Tl 

|A„| > -, Condition C2.2 holds. 

Next, in order to check Conditions C2 and C2.1, we establish some properties 

of i)k,h- 

Theorem 2.7 For any integer h, 

'4'k,h{x) = -- (^fc-l,n(a;) - ^/'fc-l,n+l(x),) (15) 
X1pk,n{x) = ^{i'k~l,n{x) +iJk-l,n+l{x),) (16) 

-^tpk,h{x) = nijjk+i,h-iix) - {h + k + l)tpk+i,hix). (17) 
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This theorem is derived directly from ( ITTj) . A recursive use of these relations results 
in the following lemma: 

Lemma 2.4 Let kQ,j,m G Z and G Z. When < ko + m — j, the function 
i'^)^'^ko,n{x) can he expressed as a linear combination of^j^o ^(x) (f = ri — m, h — 

m + 1, ... , n + m + fco — ^o) ijuhose coefficients are polynomials of n and k with degree 
not greater than m. 

Remember that the differential operator P(x, ^) is given as a linear combina- 
tion of the operators By applying Lemma [2.4[ dH]) - (ITT]) , we obtain in the 
following result: 

M 

Lemma 2.5 (Theorem 4.2 (a) (b) of [5J) Let P{x,£) = '^Pm{x){-^)"' . When 

m=0 

Pmix) {m = 0,1,... M) are polynomials and k^ < k^ — sq with sq defined as in 
([8]), {k^ G Z+, k^ G Z), the function P{x, -^)enix) belongs to V.. Then, the complex 
number 6^ := (-Be„, e^) = (P(a;, ^)e„, e^) {m,n G Z+) satisfies the following 
conditions (a) anc? (b): 

(a) : b'^ = if \m-n\>2M + ko-k^ . 

(b) : There exists a polynomial A{x) of degree not greater than M such that 

l^ml — A{n) for any m,n ^ Z"*". 

Lemmata 12.21 and 12.51 show C2 and C2.1. Thus, we have shown that the pair 
of Hilbert spaces {Lf, .(R),^^ (M)) satisfies Conditions CI, C2, C2+ and C2.1- 

C2.3. Thus, this band-diagonal matrix b"^ is illustrated by Figure [H 

Next, we point out another property of ipk,n related to Fourier series. By the 
change of variable x ^ 9 := 2 arctanx (where x = tan | ), there is an isometric map 

from the orthonormal basis system {\J\ipk,n | n G Z} of L^^j(]R) to the orthonormal 

basis system of the sinusoidal waves 

{(^e*"V e Z+} of L2((-7r, vr)). The detail 
of this relation is given in Appendix [A] The same change of variable has been used 
for a description of analytic unit quadrature signals with nonlinear phase [8] [9], 
for example. When a function passes Dini's test[l2], its Fourier series satisfies point- 
wise convergence. So, the above isometric correspondence between two basis systems 

{\J^'^k,n\n £ Z} and { ^^/^ e^"^|n G Z+} results in Theorem lA.ll in Appendix 1X1 
which shows the point-wise convergence of the expansion of any once differentiable 
function in T-i by the basis system {e„ | n G Z+}. 

2.4 Proof of Theorems [23] and [272] 

Since the 'if part of Theorem 12. II is trivial, it is sufficient to show the 'only if part 
for Theorem 12.11 Further, the the 'only if part for Theorem 12.11 is a special case of 
Theorem 12.21 Hence, we will prove only Theorem 12.21 which can be shown from the 
relation between {%) and [v) in the following theorem. 
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*: general elements 2£ 
T: non-zero elements 



vL« vL* vL* vL* vL* vL* vL* vL* 
*Y* •T* *T* 'i* 'i* *T* 'i* *T* 



Figure 1: Figure of band-diagonal matrix 6^ 

Theorem 2.8 Assume that a differential operator R{x, ^) := ^^=o ^'"(■^)(^)'" ^'^'^ 
rational coefficient functions rm{x) (m = 0, 1, . . . M). VFe denote the least common 
multiple of the denominators of rm{x) by l{x). For any X, we define a differential 
operator P{x, ^) := Ylm=oPm{x){-^)'^ - Xl{x) with Pm{x) := l{x)rM{x). Then, for 
any ko, there exists an integer kg such that the relations (i) =^ (ii) =^ {Hi) =^ (iv) 
{v) for f e q,^^^iR) hold. 

(^) ■ ^R,Lljmf = ^f- 

(a) : BpQ Lj (M) l2 (^x/ = 0. 

{iii) : r/ie f'-sequence {fn '■= (/, e„)^}^o belongs to V defined with with the 
quintuplet (P, L^^^^iR), {Jl ^l^k,,n,,J^=^, ^A^)^ ^k^n , }n=o)- 

{iv) : / e C^(R \ • rM)~^(0)) n anc/ x G R \ • rM)~^(0) safe/ies 

P(a;,f)/(,x)=0. 

(f) : / G C^^(M \ (/ • rM)"^(0)) n L^fe„)(M) and x G M \ (/ • rM)"HO) saizs/ies 
P(a;,|:)/(a;) = A/(a;). 

Proof: There exist an integer ki and a constant c such that (x^ + < c. 

Then, we choose k^ with satisfying the condition k^ < min{A;o — maxj„ (degp^ — 

m), ko + ki}. 

The property /cq < /cq + /ci yields that when a sequence /„ G L^^,^^-|(]R) satisfies 
ll/n||(fco) + - A/)/n||(,0) ^ 0, + ||P(a;,^)A||(,0) ^"o. Thus, the 

relation (i) ^ (ii) holds. 
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Since the set of singular points of the differential operator P{x, ^) is {I-tm) ^(0), 
the property < ko — niax^ (degpm — m) allows us to apply Theorem 12.31 to the 
quintuplet {P, Lf,jR), ^j^^^^^^ J^^,, L^^^,^^^^ Then, the 

relations (ii) =^ {in) =^ (iv) hold. Since the relation {iv) =^ (v) is trivial, we proved 
the desired arguments. 



2.5 Relationship to the algorithm 

The basic framework for the proof of regularity given in this paper is the same as the 
framework for the algorithm proposed in [5j and [6] which yields all the solutions in 
L^^^^^(R) of higher-order ODEs using only the four arithmetical operations on integers 
when the ODEs have no singular points or they are Fuchsian. (Note that, even if the 
value of a solution is ambiguous at the singular points of the ODE, the ambiguity 
causes no problem for discussions in a Hilbert function space where functions are 
defined by equivalent classes of the quotient space by the subspace of null functions.) 
This algorithm is based on the matrix representation of the operator -Bpa-hw* with 



respect to the basis systems {y^^fccn^jj^ {n G and {y^tp^^o ^ ^ | n G Z+} 

under the choice of spaces H = L'^f,^^{M.) and = L^^o^(]R) with k^ < k^ — sq. 

In this context, the proofs given in this paper can be interpreted as proofs of the 
validity of this algorithm, which guarantee the one-to-one correspondence between 
the square-summable vector solution of the corresponding the band-diagonal-type 
matrix-vector equation (simultaneous linear equations) and the true solutions in % of 
the corresponding differential equation, i.e., the one-to-one correspondence between 
the vectors 'm.V r\ £^(Z+) with V defined in ^ and the functions in {/ G C^'' {R \ 
S)^H\^xeW\S, P(A; X, £)f{x) = 0}. 

Since the matrix-vector equation — corresponding to the ODE is 



n 



infinite-dimensional, we should be careful of whether or not the vector corresponding 
to any solution in L^^^^(R) of the ODE P(x, ^)/ = A/ (on MyS*) always satisfies the 

matrix-vector equation ^^^^/n — 0- ^^i^ context, the proofs of {iv) =^ {ii) =^ 

n 

{Hi) can be regarded as the proof of the validity of the matrix-vector representation 
of the ODE. 

On the other hand, in infinite-dimensional case, all the solutions of the matrix- 
vector equation do not necessarily correspond to the true solutions of the ODE. 
Actually, as is shown in [5], there are vectors in V which do not correspond to any 
true solution in L^^^^(M) of the ODE (on R \ S"); nevertheless there is no such vector 
in n £^(Z"'"). However, in [5j, the statement {Hi) =^ {iv) is assumed only as 
a condition, which is C4 of [5], and its proof is omitted in that paper. In this 
context, the proof of {Hi) =^ {iv) can be regarded as a proof of the non-existence of 
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extra solutions in L^^^^(R) in our method which do not correspond to any solution 
in L^^^^(M) of the corresponding ODE (on M \ S*). In the proposed algorithm, we 
utilizes a method for the removal of the non-square-summable components from the 
vectors in V, and hence we can obtain approximations for only the true solutions 
L^^^^-|(M) of the differential equation with high accuracy. 

Thus, the proofs in this paper guarantee also the one-to-one correspondence be- 
tween the functions obtained by this integer-type algorithm and the true solutions in 
L^^^-)(M) of the differential equation (on M \ 5) . From this point of view, this paper 
contains the proofs of some propositions required in [S], which was omitted there. 
They will be given in Subsections 13.21 of this paper. 

More generally, for non-Fuchsian cases, we can use at least the statement {iii) =^ 
(iv) of Theorem l2.3[ This statement guarantees that all the solution obtained by this 
algorithm approximately coincide with true solutions of the differential equation in 
any open interval between two adjacent singular points of the ODE, and it guarantees 
that this algorithm yields at least all the solutions in C7(M \ 5*) fl {dom Bp x^-h^-ho) of 



the ODE (onR\^). 

3 A 'kind of smoothing operator' and Condition 



3.1 A 'kind of smoothing operator' for blurring endpoints 

In order to show C3, we have to check whether the contribution of the difference 
terms between two endpoints in the 'integration by parts' vanish or not as the end- 
points tend to ±oo. Usually, for functions in a Hilbert space in general, it is difficult 
to show this vanishing by a direct method because the normalizability does not al- 
ways imply smooth decays for large \x\ but may possibly allow long-lasting sparse 
oscillations with undesired peak amplitudes. For the proof based upon this vanish- 
ing, here we will introduce a convenient operator T which 'blurs' the two endpoints. 

Definition 3.1 On a space in general of locally integrable functions, define the linear 
operator T by 



Lemma 3.1 The operator T defined above satisfies the following properties: 

(Tf ) (x) is (m + l)-times continuously differentiable in M\{0} if f{x) is at least 
m-times continuously differentiable in ]R\{0}. Moreover, 



C3 





(if X = 0) . 
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(Tg) (x) = (T/) (ex) if g{x) = f{cx) (c : nonzero real constaii£^) 
lim (T/)(x) = lim f{x) if ^ lim /(x), (19) 

l(r/)(x)| < (r|/|)(x)<(T|^|)(x) (20) 

if \f{u)\ < \g{u)\ holds for \x\ < \u\ < 2|x| . 

Here we omit a discussion about differentiability at x = 0, which has nothing to 
do with the proofs in this paper. The proof of this lemma is derived directly from 
the definition of T, where the negative sign cancels out when x < because then 
X > 2x. The property (fT9|) in Lemma 13.11 is very important for our purpose because 
it results in the following lemma: 

Lemma 3.2 (In the following, denotes (^)"/ for n E Z+.) Let m E Z+. For 
functions f,gE C""(M), if there exist nonnegative integers (r = 0,l,2,...,m — 1) 
such that 

lim (T"'-(/('')^(™-'^-i)))(x) = (with (/("^^(™'"~^))(x) := /("^(x) ^("^-"-^^(x) 



/oo poo 
f{x) g^''^\x)dx and / f ^^\x) g{x)dx exist, 
■oo J — oo 

then 

/oo roo 
fix) g^"'\x) dx = i-l)"' / f^"'\x) g{x) dx . 
-oo J — oo 

Proof of Lemma 13.21 Define 

Y{x):= f{u) g^"'-\u)du and Z(x) := / f^'^\u)g{u)du. 



Then, integrating by parts (which is always applicable to integrations over a finite 
interval [— x, x] ), 

W{x) := r(x) - (-l)"Z(x) 



m— 1 



5^(-l)'-(/W(x) g^'^-'-'\x) - /«(-x) i)(-x)) . 



r=0 



Since a recursive use of f|T9|) in Lemma 13.11 results in 

lim (T"/) (x) = if ^£ G {0, 1, 2, n - 1} s.t. lim (T^) (x) = , 

x^oo x— >-±oo 

with n := maxn^, we have lim (T"'(/'-'"^(7*-"*^''^"'^^-')) (±x) = for r = 0, 1, ...,m — 1 . 

r x^oo 

Hence lim (T'^W) (x) = 0. On the other hand, 

x— >oo 

/oo 
/(x) ^('")(x)rfx 
-oo 

/oo 
f^'^\x) g{x)dx. 
-OO 
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From these facts, lim (T^W) (x) = results in the conclusion of the lemma, because 

T" is linear. ■ 
There are some other properties of T, useful for the proofs, which are summarized 
in the following lemmata: 

Lemma 3.3 Let A; G Z. For any locally integrable f in L^^^(]R), with p{x) := 
f{x) , lim (Tp) (x) = lim (T \p\ ) ix) = . 

Proof of Lemma 13.31 From the Schwartz inequality, for x ^ , 
|(TH)(x)| = — / \u'f{u)\du 



< 



1 

|x| 

1 

\x\ 



2x 



\u'^f{u)\'^ du 



1 

\x\ 



2x 



Let C ■-- 



\f{u)\^ {u^ + ifdu . (If / G Lf^^iM) , C should be finite.) Then 



2x 



u^^\f{u)\'du 



< (max(l,(4,)'^)) / 

J X 



\fiu)\'{u'+irdu 



Hence, if / G Lfj^^iR) , then \{T\p\){x)\ < 



Cmax(l, 



x^ \k 



holds for X 7^ 0. 



bmce lim 



x±oo 



for any k E Z, with ( l20l) . the proof is complete. 



Lemma 3.4 For m e Z+ , if f e C^{R) satisfies lim (T™/) (x) = , i/ien 

a;— >itoo 

hin (T™+i(7) (x) = /or gix) := x£ /(x) . 



Proof of Lemma 13.41 Since 

'.2a; 



X 



2x 



- (2x)/(2x) - x/(x) - / fiu) du = 2/(2x) - /(x) - (T/)(x) 



from Definition 13.11 (fTSl) and (fT9l) , we have 



hin (T-+i^7)(x) 



lim (T-T(?)(x) 

hin [ 2 ( T™/) (2x) - ( T-/) (x) - (T^+V) (a:)) ] = 0. 



15 



Lemma 3.5 Let f G C'^^ be a locally integrahle function satisfying lim (T |p | ) (x) = 

with p{x) := f{x) , and let g{x) satisfy the following conditions (a)-(7).- 
(a) There exists a;o > such that g{x) is at least once continuously differentiable 
for all X such that \x\ > Xq 
lim sup \g{x)\ < oo 

(7) lim sup \x£ g{x) | < 00 . 

Then, for the functions hn{x) := x^^'^g{x) {-£:)"' fix) (n E the convergence 

lim (T^'+^hn) (x) = holds. 



Proof of Lemma 13.51 The proof is by mathematical induction. 

Firstly, for the case with n = (where ho{x) = x^g{x)f{x)^, from fl20|) . the 
theorem of the lemma holds, because the conditions of the lemma guarantee that 

^Xc > and > 

s.t. < —Xc and > Xc , \g{x)\ < C i.e. |/io(a^)| < C*! a;''/(x)|. 

Next, assume that the theorem of the lemma holds for n = 0,l,2,...,n'. The 
following discussion refers only to values of x such that |x| > Xq where g{x) is 
differentiable, which creates no problem for statements about the limit as x — )■ ±00. 

From this assumption and Lemma 13.41 lim (T"'^^6„/) (x) = with bn'{x) : = 
x^ hn'{x) . Here, let 

g(x) := x-'+' ( {n' + k)g{x) + xj- g{x)] ■ (&' f{x) 



Then, since 



(M(^) = {^'''^'^'9ix)&'^'f{x)) + [xi{x-'-''gix)))-[i£f 
= hn'+i{x) + q{x), 



we obtain 



lim 

X— ^±00 



(T"'+X'+i)(a:)+(r"'+2g)(x) 



0. 



Since the trigonometric inequality and the conditions of the lemma imply that 

' dx 



lim sup |(n' + k)g{x) + x4^ g{x)\ < 00, the statement of this lemma with n = n' 



and ([20]) result in lim fT"'+^g) (x) = , and hence lim fr"'+^g) (x) = by 

X— >itoo \ / X— >itoo \ / 



( |T9l) . From these relations, lim (T"^^/i„/+i) (x) = i.e. the statement of the 

X— ^-itoo \ / 

lemma holds for n = n' + 1. ■ 
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3.2 Proof of Condition C3 

In this section, we will prove the following theorem, which shows C3: 

M 

Theorem 3.1 (Theorem 4.8 of |5J) Let P{x, £) = ^Pm{x){£)"' with p^{x) : = 

m=0 

degpm 

Pm,jX^ , o-nd let Si > So (^with sq defined in (E])). Then, for the closed extension 

j=0 

B with respect to the graph norm of the operator B defined by the action of P{x, ^) 
with domain 

D{B) = {fe n 4^,)(R) I Bf e , 

and for the closed extension C with respect to the graph norm of the operator C 
defined by 

M deg p,n 
m=0 j=0 

with domain 

D{C) = {/ G n L^,„_.,)(M)) I Cf G Ll^^m , 
the following holds: 

^/GdomB and'^neZ, {B f, ^ko-si,h) = (^f, C ^ko-si,h^ ^ 

This theorem (together with results on limits of function sequences) implies that the 
basis functions of "H^ belong to the domain of the adjoint of B under the above choices 
of function spaces and basis systems. This theorem is essential in order to show the 
statement (ii) =^ {Hi) of Theorems 12. 3[ and it guarantees that the corresponding 
number sequence {/nj^o of any true solution / in {M\pjj (0)) D {dom Ap^-}{) (or in 
C*^(M) nH for the cases where Pm{x) has no zero point) of the differential equation 
always satisfies the simultaneous linear equations '^^^b'^fn = (m G Z"*"). 
Before the proof, we establish the following preliminary lemma: 

Lemma 3.6 Let k,n ^ 7^ and j, m G Z+, and define u^n '■= max(ri + A; + 1, — ri, k + 
1). Then, for the function \^^]n{x) := {-^)"^ (^x\x'^ + l)^tjjk,n{x)^ , the function 

Rm,j,k,n{x) '■= (a;^ + 1)'''="+™- X^J^\{x) is a polynomial in x and its degree is not greater 
than 2z/fc ji + m + j + A; — 1. 

Proof of Lemma 13.61 From the definition ffTTl) of 'ipk,h{x), the function (x^ + 
^Y''"'ipk,h{x) is a polynomial in x and its degree is 2h'kh — k — 1, because the de- 
grees of the factors {x ± i) in the denominator of ipk,n{x) are not greater than Ukn 
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and the difference between the degree of the numerator and that of the denominator 
of '4'k,hix) is A; + 1. Hence, the function Tm,j,k,hix) := {-^)™'x^{x'^ + lY''*^ipk,hix) is a 
polynomial in x and its degree is 2i>kn — m + j — k — 1 when m < 2vkn + j — k — 1, 
while Tmj,k,n{x) = when m > 2z/fc + j — k — 1. 

On the other hand, the function r^,fc,n(a;) := (a;2+l)™+Vfc,ii(x) {{£r{x^ + If'''''') 
is a polynomial in x, because + l)''*""'^ contains the factor (x^ + lykh-m-k 

when m < Ukn — k. Here, when m < Ukh — k, the degree of (x^ + + 
lykfi-k (which is a polynomial) is m. When m > Vkh ~ k, the degrees of (x^ + 
l)'^^''iljk,h{x) (which is a polynomial) and (^)'^(a^^ + 1)'^'=""*^ are 2m + A; — 1 and 
St'fcri — 2k — m, respectively. From these facts, we can easily show that the degree of 
Tm,k,h{x) is 2ukn + m- k -1, 

Since Rm,j,k,h{x) = (x^ + l)'^~^"^Tmj,k,n{x) — x^Tm.k.n{x), the calculations of the 
degrees of polynomials 

2{m + k) + {2iykn - m + j - k - 1) = j + {2vkh - k - I) + m 
= j + {2vkh + m-k-l) = 2ukn + m + j + k -1 

lead us to the statement of the lemma. ■ 
By means of the lemmata in Section [3] about the operator T and the above Lemma 
13. 6[ the proof of Lemma 4.8 of the paper [5j is constructed as follows: 

Proof of Theorem O For X^J^'l.^^^^ix) := (£)"^ (^x^{x^ + if'^-'^^PkoA^)) , 

Lemma 13.61 implies that there exist finite i^', ,^ > such that | A^™^^_^^ ^(x) | < 

ir(yx2Tl)^''"'^"'"+^'~^for |x| >^ i.e. \ {x^ (x^{x^+l)''''~ '^tpko- 



K 



si , n 



< 



(Vx2 + 1 )ko+si+m-j+l 



for |x| > ^, Hence, there exists a real number K' such that 



M degpm 

m=0 j=0 
M deg p 

m 



si,n 



X 



m=0 j=0 

dog Pm. 

< E i^^'j^i 

j=0 



K' 



( + 1 fo+si+m^i+i 



for |x| > ^. 



Since Si+m— j > sq+tti — degpm > is satisfied for j < degp^, it is easily shown 
that 



M deg p. 

E E (-i)'"p;;^(a;' + i)"'"(^)™(^''(^' + i)'''"^^v^fco-..n)(x) 

m=0 j=0 



X' 



< OO 
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from the above inequality, i.e., ipko-si,h £ D{C) . Hence, Cipko-si,h is well defined 
and 

M deg pm 



m=0 i=0 

(In the following, the suffixes for j, m, ko, Si and h are often omitted if unneces- 
sary for simplicity.) 

Let f eD{B). Then, for 



Z{x) := / {Bf){u) ^ko-jA^) (m^ + rfw , 

J —X 

the convergence lim Z{x) = [Bf, V'/co-j,n holds because 5/ G H = L(^ko~si) 

and ipko-si,h e Next, define 



Y{x) := / f{u) {CiJk,-s„n){u) (x^ + l)'' du 

J —X 

M dogpm „x 



m=0 j=0 

where the convergence lim Y{x) = ( /, Cipko-j n\u ^ holds because C^/'fco-si n £ 
L^^^j(M) and / G L^^^^(]R). Then, integrating by parts (which is always applicable to 
integrals over a finite interval), 

M degpm 

Z{x) = W{x) + Y{x) with iy(x) := E p„jW^j(x) and (21) 

m=0 j=0 

m— 1 
r=0 

Here, by a recursive use of (IT^ . 

Jim (T™Z) (x) = (5/, Jim (W) (x) = (/, C^.^.,,,) (,^^.(22) 

In the following, we will show how the contribution of W{x) in (122|) behaves as 
X — )■ oo under the 'blurring' of x by the operator T defined in Section [31 From 
Lemma [3.6[ there exists a polynomial R{x) of degree not greater than 2z/fcgj^ + m + 

j + /i;o-si-r-2 such that a5.;^;:~^^^^(±x) = ^ ^)2vfepl+m~r-i where z/fc^ has also 
been defined in Lemma 4.2. Hence, with Q{x) := a;2'^feo"+™+-?'+'^0"*i"'""2 -^j^ich 
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should be a polynomial of x of order not greater than 2ukg n + + J + ~ si — r — 2, 

we have X^JT''^,%(±x) = (±xfo-"'+^-'^+^ . T for x ^ 0. Here 

note that 



lim 



'1 _|_ J_')2i^fc„ii+m-r-l 



< OO , lim 



X ■ 



(fco)V^^) -^(fco+J-m-si)( 



On the other hand, since / G D{B) C Lfko)(^) ^ ^(fco+j-m-si)(^) due to h + j - 
m—si < ko, Lemma [3]3] implies that lim (T\p\){x) = Oforp(a;) := x''°~^^~"^~'^^ f{x). 

X— >-itoo 



Then, since / G D{B) C C^-'^(M), we can apply Lemma 1331 for g{x) = ~T^_^ i \2uk ■■+m-r-i 



with kQ + j — m — Si instead of fco and with r instead of n, where p (x) = p (x) and 
/ir(x) = {^^Pko]rj\{x)^ ■ {{-^y f{^)^- (Here note that g{x) is defined for each fixed 
r, though it depends on r.) Its result 

^lirn^ [T^^\) (x) = for g(x) := (aJX-I.^U^)) " {{±rfi< 

with the definition of Wmj in (122|) implies that lim {T™'Wm j) (x) = and hence 
lim ( T"^W) (x) = . This convergence, together with the convergences (122|) re- 

X— >itoo 

suits in the required statement [Bf, ipko-j,h) i^f^^^_j-^ = (/, ^'^ko~si,n) f^j^^y because 
lim ( ( T^Z) (x) - ( T™!^) (x) - ( T^F) (x) ) = is shown from ([I9]) and ([22]). ■ 

4 Proof of {ill) =^ {iv) under CI, C2, C2+, and 
C2.1-C2.3 

In this section, we will prove that any square-summable vector / satisfying b]^fn = 
corresponds to a true solution in {R\p^{0)) fl "H of the differential equation 
P(A;x,£)/ = 0, under CI, C2, C2+ and C2.1-C2.3 . In order to show this, we 
have only to prove the following theorem and the following Theorem. 

Theorem 4.1 Assume that the quintuplet (P(A; x, {e„ }^0' "^^^ {^n}^=o) ■^o^- 

isfies Conditions CI, C2, C2+, and C2.1-C2.3. Then, any sequence f G yn£^(Z+) 
satisfies the following. There exists a function Lp G C*^(R\p^/(0)) swc/i i/iai 



P(A; X, v?(x) = and lim ^ /ne„(x) = v?(x) (23) 



n=0 



/orVxG RW(0). 
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The proof of this theorem will be constructed in this section. Theorem 14.11 implies 



N 



that /„e„ converges to a true solution of the ODE as — ?■ oo for any f G V 



n=0 



in the sense of point- wise convergence except at the zero points of Pm{.x). Thus, it 
shows that the statement {iii) =^ (iv) holds under the condition in Theorem 12.31 

Especially when pm{x) has no zero points, Theorem 14.11 guarantees the conver- 
gence with respect to the "H-norm by means of the following lemma: 



N 

Lemma 4.1 If there exists a function ip E C^^(M) such that lim y fn&nix) = ip{x 



n=0 
N 



holds for any x G M for a sequence {fn}'^=o ^ ^^(^^); then lim CS /„e„) —cp 



0. 



n=0 



n 



This is just the same as Lemma 3.10 of our preceding paper [5], and the proof is given 
in that paper. 

To prove Theorem 14. H with the projector P„ on L^^^^^(]R) to its subspace H^'^^ := 

n 

span(eo, ei, . . . e„), we will analyze the behavior of PnU = ^^l/rCr for y G yn£^(Z+) 

r=0 

as n — )■ oo. Since rj = Pnf is a solution of the inhomogeneous differential equation 
P(A;x, -^)r] = Qn with := P{X;x, ■^)Pny tautologically, we can utilize a kind of 
'continuous' correspondence between the inhomogeneous term and the solution 
rj. There, even though gn does not converge to with respect to the L^-norm, the 
convergence of ?7 to a true solution of the homogeneous equation P(A; x, ^)/ = 
can be shown with the help of the characteristic equation of A^(x, ^) in C2.2 under 
some modifications. 

Before giving the proof of Theorem l4.lt we will provide some preliminaries. First, 
in order to describe the correspondence between g^ and the 77, we will show some 
properties of the Green function for the first-order standard form of a Mth-order 
differential equation, for any intervals between adjacent zero points of pm{x), as 
follows: 

AI 

When an inhomogeneous Mth-order differential equation ''^^Pm{x){-^)"^V = 9 

with polynomials Pm{x) {m = 0, 1, M) satisfies the condition that '^x G /, Pm{x) 7^ 
with an open interval / = {z, z) and the function g{x) is continuous, we use the 
following standard form 

= M(x) f]{x) + g{x) (24) 

dx 
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with the M-dimensional vectors 



ivix)),:=^vix) (£ = 0,1,...,M-1) , (gix)), 



and the M x M-matrix 

1 



9[X) 
Pm{x) 



(if < £ < M - 2) 



mix) :-- 











Pm[x) 



Pm{x) 



P2[X) 

Pm{x) 



Pm~2{x) Pm-i{x) 
Pm{x) 



Pm[x) 



Note that [M(x)] 



for X G / if i' > i + 2. From the existence theorem, 



the m-dimensional vector- valued first-order differential equation ( 12^ has M linearly 
independent continuous solutions, because all the elements of M are bounded (hence 
Lipschitz continuity of the right hand side with respect to ff can be derived) and 
continuous with respect to x and g{x) is continuous with respect to x for x E I 
under the condition that pm{x) has no real zero. Therefore, under a choice of the 
basis vectors, there are M continuous solutions f]o{x), f]i{x), ...f]M-i{x), which satisfy 
the initial conditions (^/m(0)f = ^me (^ = 0, 1,...,M — 1; m = 0,1,...,M — 1). 
Corresponding to this, consider the following vector-valued standard form of the 
corresponding homogeneous equation P(A; x, 4^) f = : 



-^/(x) = Mix) fix) 
ax 



(25) 



Here fix) is an M-dimensional vector- valued function of x in standard form defined 
by {f{x))f, = {^xY f{^)'i it is distinct from / G £^(Z+) used in other parts of this pa- 
per. Let foix), fiix), ...fM-iix) ix G /) be its M continuous solutions which satisfy 
the initial conditions (^fmiO) = ^mi (^ G /; £ = 0, 1,...,M-1; m = 0, 1, M- 1), 
whose existence is guaranteed in a similar way to the case of (1241) . 

Define the M x M-matrix $(x, ^) by [^ix;^)]^^ := (^fmix)^ for x, ^ G / , 
which satisfies ^$(x;^) = M(x) $(x;^) and $(^;0 = for x,C, el. As is well 



22 



known, satisfies the reproducing relation 

x') ^{x'- = (x, x', ^ G /) (26) 
and another partial differential equation 

^$(x; = -$(x; M(e) (x, ^ G /). (27) 

Partial differentiability of $(x; ^) with respect to ^ is easily shown from the discussion 
about the difference under an infinitesimal change of ^, because (I27|) is derived from 
the differentiation with respect to x' of both sides of the above reproducing relation 
fl26l) and the regularity of the matrices is guaranteed by the linear independence of 
the columns. 

Here, we state a lemma about the higher-order partial derivatives with respect 
to ^ of [ $(2;; ^) ] Q especially aX ^ = x, which will play an important role later. 

Lemma 4.2 Let x,^ El. [ $(x; C,)]o m-i ^'^ partially differentiable with respect to ^ 
infinitely many times for ^ < x, where partial differentiability with respect to C, for 
^ < X includes the existence of finite partial differential coefficients from the left at 
^ = x, 

Proof of Lemma 14. 2t Since [M(x)]£^, is differentiable with respect to x infinitely 
many times for x, ^ G /, mathematical induction on m by a recursive use of (127|) 
results in the following (*) for m G Z+: 

(*) [ $(x; ] s-^^ partially differentiable by ^ for ^ < a; (x, ^ G /). 

■ 

With $(x; ^) defined above, as is well known, the relation 

^{x;x')g{x')dx' (x,^G/) 

holds with (Im)m' '■= ^mm'- Heuce, the solution fjj. of fl24|) with the initial conditions 
r/(0 = r is 

/X 
$(x; x) g{x) dx (x, ^ G /). 

Hence, if we redefine g{x) by extending its domain to M by 

( \ — \ ^ ^ ^ 

^^^^~\0 (ifxGR\/), 

under C2+, the solution of the inhomogeneous differential equation P(A; x, -^)rj = g 

for X G / with initial conditions ^^(0 = G /; i = 0, 1, ...,M — 1) can be 

written in the simple form 

r]f{x) = ($ (x; , ^) + {X^,x ,9)no (x, ^ G /) . (28) 
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with the vector x') defined by 

'$(a;; x'))^ := a;')] q , (x, x' G /; £ = 0, 1, M - 1) 

and the function 



x) 



Af-l 



(29) 



with v'^i^x) in 02"*", where lj(a;) denotes the indicator function for the interval J. 

Here, we state a prehminary lemma related to this function, where M is the order 
of P(A; X, ^) and / = (2;, i) be an open interval in which Pm{x) has no zero points. 

Lemma 4.3 Let ^ E I . Under C2+, C2.2, and C2.3, for any u E I greater than 

^, ^K^ u > and ^ric G Z+ such that | ( > 6n)->^o | — — If -^^^ '^'^^ n G Z+ 
greater than ric- 

Lemma 4.4 Under C2.2, [N{x, ^))*^ can be expressed as the finite sum 

M 



m=0 

where the functions z/^ (m = 0, 1, . . . , M) belong to C^{R). 

The proof of this lemma follows easily by mathematical induction on M. 

Lemma 4.5 Under C2.2 and C2.3, for any real numbers a and b in I such that 
a < b, a function f G C^\l) satisfies the relation 



^CabER S.t. G Z+, 



<Ca> 



Proof of Lemma 14.51 Under C2.2, (^N(x, -^))^^ is well defined. From Lemma 
m with := uUx)f{x) G C^'\i) C C^'la,b], 



M „6 



r fix) {{N{x,±)f'el{x))dx = Y, I /(^) (^^„.(a:)(e^)(™)(x))dx 
J] (-ir / li^;:\x)ei{x)dx 

m— 1 



m,=0 



r=0 
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Since /i^ ("^ = 0, 1, ... M - 1) and (n G Z+) belong to C*^[a, 6], the functions 
fim{x) and em''(a;) belong to C^[a,b] for r < M. Hence, under C2.2 and C2.3, 
all the maxima := max and A^"^^ := max |a*-^-'(x)| {0 < m < M and 

< r < M) are finite, with a{x) in C2.3 . From these facts, 

M / m-l 



fix) [{N{x, ^)Y'ei{x)) dx\<J2[ib- a)M(:")i(°) + J] 2M«i( 

m=0 \ r=0 

where the right hand side is finite and does not depend on n. 



such that f^,u{x) 



31 Let ^, M G / = (z, z) and let /j,^ be a function in C {!) 



^.31 


Let ^, M € 


$( 










A/-1 



for X G [^,u\. The existence of is obvious 



from the extension of the function to the intervals {z, ^) and (m, z) by the Taylor 









Af-l 



expansions of ^=^= up to the M-th order term about x = ^ and x = u, 

Pm{x) 

respectively, because of Lemma 14. 2[ 

Under C2+, C2.2, and C2.3, since N{x, e^{x) = A„e2(x), 



Af(xc,«,e<>>^, = Af / fU^)el{x)dx= j f{x)({N{x,±)fel{x))dx 
This and Lemma 1^31 result in ^C^,^ G M such that |A„|^ |( , 6n)->^o| — The 

I '^'^ I 

condition liminf > in C2.2 implies that there exist an integer ric and a positive 

n— >oo n 

constant c such that |Ari| > cn may be guaranteed for any n greater than ric- Hence, 



C C 
I ( Xii,u , e^)^o I < j^^^ for any n G Z+ greater than n^. With := the 

lemma holds. ■ 
Next, as another tool for the proof of the theorem, we will consider the problem 
of finding the solution of the differential equation P(A; x, -£:)r] = g for x E I = {z, z) 
under the constraints rj{xj) = tj (j = 0, 1, M — 1) for a sequence xq < xi < ... < 
xm-1 in I, instead of giving the M initial conditions only at x = ^ {z < < xq). 
For this problem, define the M x M-matrix T by (T)^.^ := fm{xj) {j = 0, 1, M — 
1; m = 0, 1, M — 1) with the solutions fm ijn = 0, 1, M— 1) of the homogeneous 
differential equation P{\; x, f = for x G / = {z,z) with the initial conditions 
/(^) = 5trn = 0, 1, M — 1) where z < ^ < xq. Then the following lemma 
holds concerning the invertibility of T: 

M 

Lemma 4.6 When P(x, ^) = '^^Pm{x){-£:)"^ for polynomials Pm {m = 0, 1,...,M) 

m=0 

satisfying ( x G /, Pm{x) ^ O), for any y E I not smaller than ^, there exists a 
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sequence of finite intervals [ao,&o], [^i^^i], • • • [oa/-i, ^Af-i] with y < ao < bo < ai < 
hi < ... < ctAf-i < &Af-i < z such that T may he invertihle when Xj G [aj,bj] 
(j = 0,l,2,...,M-l). 

Proof of Lemma W^ Define nxn-submatrices T„(xo, Xi, x^.i) {n = 1,2, M) 
by (f „(xo,Xi, ...,Xn^i)).^ := fmixj) (j = 0, 1, ...,n- 1; m = 0, 1, ...,n- 1). Then 

Tm{xo,xi, ...,xm-i) = T. Since the statement that fo{x) = for any x in {y,z) is 
contradictory to the uniqueness theorem and the initial condition at x = there 
exists Xq such that y < Xq < z and /o(a;o) 7^ 0. Then detTi(a;o) = fo{xo) 7^ 0. 
From this initial statement, we can carry out the following mathematical induc- 
tion: When detTj{xo,Xi, ...,Xj_i) 7^ 0, there should exist Xj in {xj_i,z) such that 

det Tj_|_i(xo, Xi, Xj) 7^ 0, because det Tj+i(xo, Xi, x) = for any x in 

j 

(xj^i, z) would imply Cm{xo, Xj_i) fm.{x) = for x > Xj_i with Cj{xo, Xj_i) = 

m=0 

det Tj(xo, xi, Xj_i) 7^ which is contradictory to the uniqueness theorem and the 
initial condition at x = From this mathematical induction, there exists a sequence 
y < Xq < xi < ... < xm-1 < z such that det TAf(xo, xi, xm-i) 7^ i.e. det T 7^ 0. 

Next, from the conditions for P(x, ^) and the existence theorem, det T = det TMixo, xi, xm-i) 
is M-times continuously partially differentiable with respect to Xj {j = 0, 1, M — 1) 
in I and moreover totally differentiable J*^, and hence it is locally Lipschitz continu- 
ous there. Therefore, with the conventional vector notation x G M^^ defined by (x)j = 
Xj (j = 0, 1, M — 1), if det T„(xo, xi, x„-2, x) 7^ 0, there exists a neighborhood 
[/e(x) = {u \ IIm — x|| < e} (e > 0) in J*^ such that det T„(mo, mi, 7^ for 

any u G Ue{x). Since {-u | G [xj — 6j, xj + 6j] (j = 0, 1, M — 1)} C Ue{x) holds 
at least for < 6j < (j = 0, 1, M — 1), the lemma holds with aj := Xj — 6j 
and bj := xj + 6j (where z < bj < Oj+i < z is satisfied for an appropriate choice of 
sufficiently small 6j and ■ 

Under the existence of a sequence with invertible T guaranteed by this lemma, 
we have another lemma with the definition of the vector bg defined by 

{bg)j := {Xi,., , 9)n^ U = 0, 1, M - 1). (30) 



Lemma 4.7 When the sequence z < xq < Xi < ... < xm-i < z is chosen so that T 
is invertible, the solution of the inhomogeneous differential equation P(A; x, -^)rj = g 
for X G / under the constraints Tj{xj) = tj (j = 0, 1, M — 1) 
[where z < ^ < xq < xi < .... < xm-i < z) is 

%-i(r-fe,)(^) = ('^(^; 0> T-i(r- bg) ) + {xi,x , c/)ho 
with the vector te M*^ defined by (i) . = tj {j = 0, 1, M - 1). 
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Proof of Lemma \A.l\ Since the homogeneous differential equation P(A; x, ^)/ = 
for X G / is a special case of the inhomogeneous differential equation with (7 = 0, 
from f l28|) . the solution of the homogeneous differential equation P(A; x, ^)/ = for 
X E I with the initial conditions ^ /(O = (^)^ = 0,1,...,M — 1) is ff{xj) = 
r) . As special cases, we have fmixj) = 1^) with the vector 



0,1, 



Im defined by ( 1^ 
we have 



tr 



e 

M 



. Define the M- dimensional vector such that (^t)^ 



1). Since rj^ixj) — [bg] = frixj), from the above relations. 



/m(a;o) 
fmixi) 



M-l 

E 

m=0 



M-l 

E 

m=0 

^ fmi^M-l] 

Hence, we can show that the function 



Tf. 



it), 



is the solution of P{X; x, -^)ri = g for x G / satisfying the constraints ri{xj 
(j = 0, 1, M — 1) for the sequence z < ^ < xq < Xi < .... < xm~i < xm < z, where 
the uniqueness of the solution satisfying these constraints has been shown also. ■ 

By using these preliminaries, now we are able to construct the proof of Theorem 
14.11 as follows; 

Proof of Theorem 14.11 

Suppose that y G Vnf{Z+) has no function (p G C^^^(M\p^/(0)) satisfying ([23]). 



Then, the basis syetem {e^ | n G Z+} is a CONS of = L 



the function 



y := ^^^nCn belongs to Ti. With the projector P„ on "H to the subspace Ti^"^^ : 

ra=0 

span(eo,ei, ...,..„;, 



holds. Hence, there exists 



e„) , the convergence lim \\Pny — y\\ /^,o^ 

a subsequence {nu}'^=n such that lim (Pn^yMx) = y{x) (a.e.). 

Therefore, from the assumption that y G (y^^\V) C V and Lemma [4.6[ without 
loss of generality, we can show the existence of a sequence ^ < xq < xi < .... < 
xm-i < xm (where M is the order of the differential equation P(A;x, ^)/ = O) 
satisfying the following conditions (a)- (d): 
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(a) For j = 0, 1, M, the hmits hm (P„^ y){,Xj) exist and hm (P„^ y){xj) = tj. 

(b) The M x M- matrix T is invertible under the definition by (T)jm = fm{xj) 
(j = 0, 1, M — 1; m = 0, 1, M — 1) for the continuous solutions fm of 
P(A; X, £)/ = with the initial conditions {£Y /(O = Se^ (£ = 0, 1, M - 

!)• 

(c) /(xm) 7^ for the true continuous solution /(x) of P(A;x,^)/ = which 
satisfies f{xj) = tj (j = 0, 1, M — 1). 

(d) Pm{x) has no zero point in [^,Xm]- 

Define 5'„(x) := (P(A; x, ^)(P„ y)) (x) . Then gn{x) belongs to C°°(]R) because 
C2 implies that it can be expressed as a finite sum of the basis functions e^{x) 
of li. Moreover, define the M-dimensional vectors t and by := tj and 

{t'^)j '■= {Pnyy){xj) (j = 0,1,...,M — 1), respectively. Since the condition (d) 

requires that ^,Xo,Xi, . . . ,Xm belong to an open interval I in which Pm{x) has no 
zero points, from the definition, P„^ y{x) is just the solution of the inhomogeneous 
differential equation P(A; x, ■£:)f]{^) = dn^x) for x G / under the constraints ri{xj) = 
(t . (j = 0, 1, M — 1) . Therefore, from Lemma [4.71 



no 



where the function x^,xm the vector bg have been defined in ( l29l) and ([30] 
respectively. On the other hand, with g = in the same lemma, similarly we have 



/(xm)= ($(xm;0, T-if 



Hence, 



{Pn. y) (xm) - if) (xm) = ( ^xm; 0, ( (f^ - f) - bgj ) + (xc,.^^ , grJn<>- 

From the definitions, the limits bg„^^ — )■ and , S'n,^) — > as u ^ oo holds 

if the convergence lim (xf i , dn) = holds for j = 0, 1, M . 

Now, we will prove this convergence, as follows: From C2, when n > 2i, it 

n+£ i 

is easily shown that gn{x) = e^)^o e^(x) , because ^m2/r- = 

r=n— £+1 r=n~£ 

max(7i, m+£) 

(m G Z+) holds for y G fl £^(Z+) and hence e^)^o = 6^ y,. vanishes 



when m + i < n. Hence, from Lemma [4. 3[ when n > £ + 1, 

n+i 



r=n— £+1 
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r=n-e+l 



n+£ \ / A\ I 



Here, for n — i + l<r<n + 



{9n, e^)^o 



< 



n—r 

— V I 



e=-e 



From C2.1, the finite supremum K' := sup 



6" I 

I 



rez+\{o}, nez+ r 



M 



I 2/r+£ I ■ (31) 

exists. Hence, from f l3T|) . 



for n — £+l<r<'ri + £, we have 

n—r 



wo 



IVr+e 



(32) 



where the last inequahty is derived from the Schwartz inequality. From the inequal- 
ities (15T]) and (152]) . for n > 2 i, we have the inequality 



\ 



r=n-2^+l 



Since ^ is square-summable, IVr+e |^ — as n ^ oo. Therefore, 

r=n-2£+l 

lim (xco;, , 5'n)wo = and hence lim , gn,)n<' = 0- 



(33) 



Thus, we have proved that lim (xf , 9nS)u'> = > lim = and 

lim {Xi,XM ' 9n^)no = 0- These convergences, together with the convergence lim t = 
t which is identical to (a), lead us to the conclusion that lim (Pn^yMxM) = /{xm), 

which is contradictory to (c). Therefore, the assumption that y n£^(Z"^) does 
not satisfies the conclusion of Theorem 14.11 That is, we obtain Theorem 14.11 



5 Discussion 

The m-th order derivatives (m = 1, 2, . . . , M — 1) of the eigenfunction / in -Z^^^qJ 
do not always belong to L^j^^^(M). For example, for the differential operator. 



P(a;, A) = (3^2 ^ 1)2 (A) 2 + 6x(3x2 + 1) (^) - (3^^ + 1)^ - ISx^ , 

1 



an eigenfunction associated with the eigenvalue —6 is f{x) = — - cos(x^ + x) , 

which belongs to L2(M). However, -;-f(x) = — sin(x^ + x) — - — -77 cos (x^ + 

ax [3x^ + 1)'^ 
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x) ^ L^(R). There are many similar examples. In order to discuss regularity in 
the framework based on the Sobolev space, some transformation or some change of 
variable is necessary, for these cases. 

Our proof does not require any assumption about whether or not the m-th order 
derivatives (m = 1, 2, . . . , M — 1) of the eigenfunction / belong to L^^^^-|(R). Hence, 
it can be used to show regularity for these cases without any transformation or 
modification. 

6 Conclusion 

When an operator ^/j.lz^ ^(m) is defined on C^(]R) fl L^^^^(M) as the action of an 

M-th order differential operator ^) with rational coefficient functions, we have 
proved the regularity of the eigenfunctions of the closed extension on L^^^-j(M) of 

the operator ^_r,l2^ ^{m) except at the singular points of the ODE R{x, ^)f = A/, 
This derivation does not require any assumptions about the m-th order derivatives 
(m = 1, 2, . . . , M — 1) of the eigenfunction. In particular, for fco = 0, we have proved 
it for the usual L'^iW). 

The proof has been constructed in two steps: The first step has proved regularity 
in a general framework under several assumptions. The second step has shown that 
the above mentioned operator satisfies these required assumptions. 

In the first step, the differential operator is treated as an operator from a dense 
subset of a Hilbert space "H to another Hilbert space "H^ which contains "H (in the 
sense of sets), and this operator can be represented in matrix form with respect 
to appropriate basis systems of 1-i and l-C^ . The proof in this framework has been 
based on the implications (z) =^ (ii) =^ {Hi) =^ (zf ), in Theorem 12 . 3 1 with a more 
general framework: (z) the kernel of the closed extension of the operator on Ti defined 
as the action of the differential operator, (ii) the kernel of the closed extension 
of the operator from a dense subspace of Ti to "H^ defined as the action of the 
differential operator, [iii) the space of square-summable number sequences satisfying 
the simultaneous linear equations corresponding to the matrix representation of one 
of the above two operators and {iv) the space of 'regular' solutions of the differential 
equation which are continuously differentiable M times at any points except the zero 
points of the coefficient function of the highest order. This general framework was 
used also for an integer-type algorithm for solving higher order homogeneous linear 
ordinary equations in our preceding paper [5]. 

In the second step, we have shown that the choices "H^ = ^ (R) and the basis 

function systems in ([9]) satisfy the conditions required for the framework in the first 
step. 

The proofs in the two steps have easily been constructed except for two points; 
one is the proof of {iii) =^ (iv) in the first step and the other is the proof of the 
fact that the choices satisfy condition C3 in the second step. For the latter point, we 
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have developed a kind of smoothing operator, as a tool. Our proof of {in) =^ (iv) 
has been made by means of a modified kind of continuity of the solutions of an 
inhomogeneous equation with respect to the inhomogeneous term. 

Similar proofs of regularity may be possible even for other choices of function 
spaces and basis systems satisfying the conditions in this paper or similar type of 
conditions, which will be a topic for future research. 
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A Relationship with Fourier series 

The basis systems used in our methods are closely related to Fourier series ex- 
pansions of functions defined on the interval [— vr, vr], by the change of variable 
9 = 2 arctan x ( or x = tan | ). By this change. 



1 2 , . . .9 .... 

(34) 



dx = - sec^ - d9, d9 = — dx , + 1 = sec^ ^ , arg(x ±i) = =F- (9 — tc 



x±i = ±x e+ 2 sec - , = - e 

2 ' x + z 



and then 



Here, define 



^.o,r.(tan^) = z'=+i(-l)^ e^('^+^)^ cos'^^^ ^ . 



i(fe + l) / 




+7r) 



sec''+^ 



9 



Q 

/(tan -) ( if - vr < < tt) 

(35) 



2' 



( if ^ = ±7r) . 

Then, from the relations (l34|l . we have a kind of isometric relation 



(/, 9\k,) = / fm UO) d9 . (36) 
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The relation (134|1 and the definitions (fTTl) (135!) result in 

e™'' f if - TT < X < vr 



^k,m = <( ^ , (37) 

( if x = ±7r) 

where we have the characteristic equation 

^^k,m = ^ni>k,m i-n<e< vr) (38) 



which corresponds to the characteristic equation flMj) in the x-coordinate in Section 
[21 From ( 157]) . the expansion of / G L(fc„)(]R) with respect to the biorthonormal basis 

system {J^ tljk,h \n eZ}, 



/(x) = ^ 'fn^k,n{x) with fn=^ f (x) i!k,n {x^ + i f dx (39) 

just corresponds to 

/>) = V Fiie^"^ with h:=^ I ^ f{9)e''''U9 , (40) 

by the change of variable x ^ 9 and the relation 

(-1)" ■■ 

= S^/r. • (41) 

The correspondence introduced above provides us with the following theorem: 

N 

Theorem A.l // a number sequence {/n}^o ^ ^^(Z^) satisfies lim ( > /nCn ) — 

/ 



n=0 

= with f E C^(K), then lim > fnGnix) = f{x) holds for any x G M. 

(fco) A^^oo ^ 

n=U 



Proof of Theorem \A.1\ 

From the correspondence ([9]) between the unilateral orthonormal basis system 

{e„ I n G Z+} and the bilateral orthonormal basis system {^J^^'ko.h \ n E 1,} of H, it 

is easily shown that the coefficients fa {n E Z) in the expansion (1391) correspond to 

oo 

the coefficients /„ (n G Z"*") in the expansion f{x) = fn^nix) by the relation /„ = 

n=0 

ffik^^- Hence, the condition / G L^^^^(M) is equivalent to the conditions {/„}^o ^ 
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£2(Z+), {/ii}^=„oo e £2(Z) and {{-lfFn}^^_^ G ^^(z), where the equivalence to 
the last condition can be shown by (14T|) . Hence, under the condition / G L^^,^^^(]R), 

the function /(^) defined in ([35]) for /(x) in this lemma belongs to 7r,7r]) C 

f^"" \f(u)\ [^'^ \hu)\ 

TT, tt]). This fact implies that / c/m < oo and / du < oo if 

J±^T^ d-u ^ e + u 

< e < I and |^| < vr — 2e. Moreover, the function f{6) is continuously differentiable 
once in the interval [— tt + e, tt — e], from fl35|) and / G C^(M). 

From these facts, for ({>e{t) '■= f{6 + + f{6 — t) — 2f{9), it is easily shown that 



/ dt < oo for any 6 G [— tt + 2e, tt — 2e]. Thus, Dini's test is passed for 

Jo ^ 

point-wise convergence of the Fourier series, for any 6 G [— 7r + 2e, vr — 2e]. Hence, from 

N N 

dSTD and (gOD, lim ^2 V (-1)"^^^ ii(^) = lim V e*"^ = f{e) holds for 



h=-N h=-N 

-N-1 

any 9 G (-7r+2e, 7r-2e) for any < e < f . Moreover, lim a/2 Fn'ijjko^r. 

ri=— TV— fco— 

(rf = 1, 2) for any 6 E {-ir + 2e, vr - 2e), because {(-l)"Pk}^^_^ G (as is 

mentioned above) and 

n=—N—ko — d h=—N—ko—d h=—N—ko—d 

1 r -N-1 

< ^ii±A_ ( j2 m 



(1 + x2)^'0+l 



TV 

Therefore, lim V2 V ii(^) = f{e) {d = 1,2) for any 6 e {-n + 

ri=— A''— fco"*^ 

N 

2e, vr - 2e). This fact and (|41]) imply that lim ^feo,n(^) = /(^) = 1' 2) 



Af-s>oo 

h=—N—ka — d 



, -7T + 2e vr-2e\, ^ ^ ^. ±vr ^ 2e 
for any x G tan , tan tor any < e < t- oince lim tan 

' 2 2/ 2e^0+ 2 



N 

±oo, lim > fh'ipkoh{x) = f{x) {d = 1,2) holds for any x G M. Since the 

n=—N—ko—d 

N 2N+ko+d-l 

'matching' in ([9]) results in fh'4'ko,n= /nCn (ci = 1, 2) where the 

h=—N—ko—d n=0 

last equality should hold because {e„ | n G Z+} is a basis system of Ti, the conver- 

2Af 2Af+l 

gences lim 7^/ne„(x) = f{x) and lim fn^nix) = f{x) hold for any x G M, 



n=0 n=0 
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N 

and hence lim fn^ni^) — fix) holds for any x e R. ■ 

AT^oo — ' 
n=0 
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